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■irAisfiuef^  ~  - — - 

For  miiy  cnoisioiied  ippBcatioiii  of  wlrotcn  mior  nilworla  (WSNf),  the  iRfrrnetion  procoiviBg  iavolvcs  tltillng  with  distributed  d«t«  io 
the  context  of  iceunte  rignil  dctcdioB  end  energy-offlctat  routing,  whkb  bm  been  letko  rosenrcb  topics  to*  mipy  years,  respectively.  In 

this  pepe^  we  relMetbesetnw  aspects  riiJoiBtop«iaiiatiea«  Coilsidering  the  scenario  ofuslogdNttibutedrJidti^lioimori  to  detect  the 

pretence  of  sa  object  througb  active  flcasiDgR  we  fermuMta  the  proUem  of  enenfH^lont  routiag  for  sigita]  detwctiai  wider  the  IHcymaiK 
Pearson  criterion,  apparent^  for  the  first  tfane.  The  joint  optimizatioD  of  detection  and  raatfag  to  carried  oat  la  a  fhsioii  center  which 
procoBipnttt  the  rontos  as  a  function  of  the  geographic  location  to  be  moahored.  Aceordiai^,  wo  propose  three  diflhrcnt  routtog  wctrics  that 
ain  at  an  appropriate  tradeoff  between  the  detoctioa  pertiM-naace  aad  the  energy  expeiiditiire.  In  particular,  each  metric  rdates  the 
dete^on  perfsrataace  explicitly  In  terms  of  prohabUlttos  of  detection  end  MIse  alenn,  wBh  tiie  energy  couaned  in  sensing  and  routiag.  The 
routing  prohimi  are  formulated  as  confalnatorial  optimization  programs,  and  we  provide  soluthHis  drawing  oa  operathma  research.  We 
IHrescat  extensive  simulation  rcsuhi  that  dcaMinstrato  the  eaegy  and  detection  performance  tradtofEi  for  each  proposed  routing  metric. 
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Energy-Efficient  Routing  for  Signal  Detection  in 
Wireless  Sensor  Networks 

Yang  Yang,  Student  Member,  IEEE,  Rick  S»  Blum,  Fellow,  IEEE,  and  Brian  M,  Sadler,  Fellow,  IEEE 


A^sirac^For  m&ny  covlsloiied  appUcadom  of  wireless  sensor 
networks  (WSNs)^  the  informatioii  processing  Involves  dealing 
with  distribyted  data  in  the  context  of  accurate  signal  de« 
tection  and  energy-efficient  routing^  wbldi  have  been  active 
research  topics  for  many  years,  respectively.  In  this  paper,  we 
relate  these  two  aspects  via  joint  optimhoitloii.  Consideiiiig  the 
scenario  of  using  distributed  radar-like  sensors  to  detect  the 
presence  of  an  object  through  active  seosliig^  we  formulate  the 
problem  of  energy-effident  rooting  for  signal  detection  under 
the  Neyman-Feaiaon  criteiiim,  af^arently  for  the  tliiie.  The 
joint  optimization  of  detection  and  routing  b  carried  out  In  a 
fusion  Center  which  precomputes  the  routes  as  a  functkm  the 
geographic  loeation  to  be  monitored.  Accordingly,  we  propose 
three  different  routing  mietrks  that  aim  at  an  i^pn^riate  tradeoff 
between  the  detection  perloriiuiiuc  and  the  eneigy 
particular,  nietrlc  relates  the  detection  peiformance  aqdidlly 

in  terms  of  probabilities  of  ^tectk^n  ami  fiilie  akrni,  wtlh  the 
energy  consumedi  in  sensing  aidd  routia^.  The  routing  problems 
are  Ibrnnilided  as  cmnblnatoiial  oirtlmlzation  prograins,  and  we 
provide  solutions  drawing  m  operalhms  ramrciL  pin^nt 
extenrive  rimulathm  results  that  demonstrate  the  enargy  and 
^tetection  performaiH^e  tnuleol&  for  proposed  rmithigmetirk. 

Index  Terms — QuiiMnatorbl  optimization,  dioitait 

fmtii,  ei^fgy  Lagran^an  relaiiatlon,  NeymaiiHPeari^ 

critemn,  paranastric  i^rieslpath,  touting,  detection,  wire¬ 

less  semor  nriworks. 


1,  iNTOODUCnON 

WIRELESS  sensor  netwcffks  (WSNs)  have  numerous  ap“ 
plications  in  both  the  military  and  civil  domains  (see 
[1]  and  [2]  for  tumrki  treatments),  and  signal  detection  (c,g,, 
object  or  event  de^don)  often  serves  as  the  initial  objective. 
In  fact,  si^al  detection  in  WSNs  has  been  extensively  studied 
fOT  many  years  (see  [3>“[8]  and  references  therein).  The  tradi¬ 
tional  techniques  of  signal  detection  in  WSNs  or  multisensor 
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systems  usually  feature  a  fusioit-ceiiUric  appit>ach  by  assuming 
that  there  Is  a  desipiabed  fusiem  center  where  a  final  decisicHfi 
will  be  made,  and  largely  focus  on  various  deiecticHi-bMed  pet- 
formance  measures  such  as  error  probabilities,  The  limited  en- 
ergy  supply  is  a  critical  issue  that  poi^  many  challanges  to 
network  design  and  management.  Thus^  recent  research  effcsis 
have  tried  to  opdinizie  the  signal  cfetecticHi  pcrformaiK^  jointly 
with  the  energy  dissipation  in  the  systom.  Ibr  example,  [9]  pro¬ 
poses  a  sensor  censoring  schen^  for  distriboted  detoction  to  re¬ 
duce  the  colt  of  transmission,  where  each  sensca'  transmits  only 
highly  Informative  chservations.  Along  this  line,  [10]  and  [II  ] 
further  propose  different  but  effective  measurement  car  transnus- 
rion  approaches  to  optimii^  the  cost  and  detection  pertormance. 
A  hybrid  energy-driven  detection  scheme  is  proposed  m  [12] 
where  the  transmissIcHn  from  each  sensor  may  vary  from  Lbit 
decision  to  all  its  observations.  Three  options  for  multihop  fu- 
siem  are  proposed  in  [1 3].  It  is  noticed  that  [9]-[  1 1 J  all  assume  a 
simple  pardlel  sensor  network  topology,  which  however  could 
reiailt  in  significant  energy  expcndituie,  especially  when  some 
sensornodbsarefir  firomthefusioncenter.  Incontrast,  both  [12] 
and  [13J  consider  multihop  information  transmission,  although 
multihop  routing  is  not  axitessed. 

The  constrained  energy  supply  in  typical  WSNs  further  ren¬ 
ders  the  task  of  finding  and  maintaining  energy-efficient  routes 
nontrivial*  (Consequently,  it  is  highly  desirable  to  find  methods 
for  energy-efficient  route  discovery  and  relaying  of  data  from 
the  sensor  nodes  to  the  fusion  center  so  that  the  lifetime  of  the 
network  is  rnaximized*  Over  the  past  few  years,  many  new  al¬ 
gorithms  have  been  proposed  for  the  routing  problem  in  WSNs 
(see  [14]  and  references  therein).  For  example,  consider  the  ap¬ 
proaches  in  [15]  and  [16],  The  authors  of  [15]  propose  a  min¬ 
imum  cost  path  routing  algorithm  which  aims  at  maximizing  the 
networic  lifetime.  Alternatively,  the  authors  of  [16]  introduce  in¬ 
formation-directed  routing  with  the  objective  to  minimize  com¬ 
munication  cost  while  maximizing  information  gain. 

Obviously,  the  literature  on  energy-efficient  routing  or  signal 
detection  in  WSNs  is  abundant,  respectively;  however,  there  is 
much  less  research  on  routing  for  signal  detection,  let  alone  the 
consideration  of  their  joint  optimization*  Some  recent  related 
work  includes  [17]  and  [18].  In  particular,  [17]  studies  a  se¬ 
rial  fusion  method  for  collaborative  signal  detection  in  WSNs, 
and  also  develops  a  space-filling  curve-based  routing  technique, 
where  an  imaginary  curve  is  drawn  in  the  geographical  region 
of  interest  that  fills  that  region,  and  the  sensor  network  is  passed 
across  along  this  curve  to  implement  serial  fusion.  In  [18],  the 
authors  consider  the  problem  of  routing  for  detection  of  a  corre¬ 
lated  random  signal  field,  and  further  propose  a  Cherrofifiouting 
metric  using  the  Chemoffi  information.  This  routing  metric  ne¬ 
cessitates  further  consideration  of  practical  is^^sui^,  nevertheless 
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it  is  very  interesting  and  attractive,  since  it  allows  one  to  char- 
actoize  the  overall  detection  performance  based  on  the  sensojrs 
along  a  route.  In  contrast  to  both  [17]  and  [  1 8],  we  considbr  here 
a  much  different  signal  detection  and  routing  problem  where  an 
emitter  or  reflector  produces  a  signal  of  interest  from  a  point  in 
space,  and  we  wish  to  find  optimal  routing  which  strives  for  a 
balance  between  the  energy  consumption  and  detection  perfor¬ 
mance.  More  specifically,  we  focus  on  energy-efficient  rcRiting 
while  optimizing  the  Neyman-Pearson  detection  critericm. 

The  Neyman-Pearson  criterion  [19),  as  an  altemadve  to  the 
common  Bayesian  framework,  is  fundamental  to  the  develop¬ 
ment  of  the  theory  of  hypothesis  testing,  and  It  is  commonly 
employed  for  surveillance  systems  [20],  [21],  Recognizing  a 
basic  asymmetry  in  the  importance  of  the  hypotheses  (i.e.,  the 
target  present  hypothesis  Hi  and  the  null  hypothesis  the 
Neyman-Pearson  criterion  maximizes  the  probability  of  detec¬ 
tion  subject  to  a  predetermined  bound  on  the  probability  of  false 
alarm.  Thus,  it  will  be  advantageous  to  employ  this  criterion 
in  WSNs  tasked  with  solving  signal  detection  problems.  How¬ 
ever,  to  the  best  of  our  knowledge,  th^  has  not  been  ai^  work 
on  energy-effident  routing  for  detection  that  specifically  con¬ 
siders  the  N^man-Pearson  entmon.  This  basically  motivates 
our  work. 

In  this  paper,  we  address  the  problem  of  energy-efficient 
routing  for  signal  detection  under  the  Neyman-Pearson  crite¬ 
rion.  We  consid^  a  set  of  sensors  which  are  distributed  over 
an  area  to  delect  the  possible  presence  of  an  illuminated  taiget 
at  a  specific  location  through  active  sensing.  The  ftision  center 
is  assumed  to  have  the  knowledge  of  geogr^Ucal  location  of 
each  sensor,  and  will  determine  the  location  to  be  probed  as 
well  as  the  communications  routing  (i.e.,  centralized  routing  is 
employed).  The  observations  at  each  sensor  are  assumed  to  be 
independent  conditioned  on  the  hypothesis.  Each  sensor  node 
makes  a  measurement  to  test  for  a  target  at  the  predetermined 
location,  and  then  transforms  the  measurement  into  a  likelihood 
r^o.  This  can  be  regarded  as  one  particular  type  of  in-netwoik 
processing  in  WSNs;  diat  is,  instead  of  sending  the  raw  data, 
each  sensor  node  uses  its  processing  ability  to  transform  its 
measuremimt  into  a  sensor  likeiil»>od  ratio,  a  sufficient  statistic 
for  the  problems  we  consider  (conditionally  indqiendent  sensor 
data)  that  provides  significant  compression.  After  collecting  all 
the  likelihood  ratios  from  the  sensor  nodes  along  a  route,  the 
fusion  center  will  make  a  final  decision. 

Our  work  exploits  a  synergistic  approach  by  incoiporMing  as¬ 
pects  of  signal  processing,  computer  science,  and  operations  re¬ 
search,  and  can  be  regarded  as  om  special  type  of  the  cross-layer 
network  design.  To  summarize,  the  novelty  and  contributions  of 
this  paper  are  threefold. 

1)  We  formulaK^  the  problem  of  enogy  efficient  routing  for 
detection  under  the  Neyman-Pearson  criterion,  widely 
adopted  fer  target  detection  and  surveillance  related  ap¬ 
plications.  This  formulation,  as  far  as  we  are  aware,  is  the 
first  one  which  accounts  for  both  the  energy  consumption 
in  sensing  and  routing,  and  delation  performance  (in 
terms  of  detection  probability  and  false  alarm  probaWlity) 
at  the  same  time. 

2)  We  propose  three  routing  metrics,  and  for  each  of  them, 
the  ^tection  performance  and  the  energy  expenditure  are 
cof^idered  jointly  in  a  different  but  interesting  way  by 


imi 


which  an  aypropriate  tradeoff  between  them  1$  attained.  We 
also  ftmnulaie  the  routing  problems  under  th€^  metrics 
into  different  combmatoria!  optimization  programs,  where 
the  objectives  are  to  maximize  or  minimize  a  quantity  of 
interest. 

3)  We  provide  algorithms  for  solving  those  formulated  integer 
programming  problems,  based  on  state-of-the-art  opera¬ 
tions  le^Mch  results.  In  cases  where  the  optimum  solution 
is  computationally  intensive,  we  provide  methods  where  * 
we  can  specify  the  acceptable  complexity  and  a  solution 
that  can  be  made  arbitrarily  close  to  the  optimum  solution 
when  its  complexity  is  increased  appropriately.  We  also 
present  exteaisive  simulation  results  which  help  quantify 
the  difference  between  these  proposed  routing  meirics. 

The  remainder  of  this  paj«r  is  organized  as  follows.  In  Sec¬ 
tion  E,  we  formulate  the  problem  of  routing  for  signal  detec¬ 
tion,  which  includes  the  description  of  the  system  and  signal 
model,  the  analysis  of  both  Neyman-Pearson  detection  and  en¬ 
ergy  consumption  in  WSNs,  and  the  introduction  of  a  graph 
model  In  Section  III,  we  propose  the  maximum  path  detec¬ 
tion  efficiency  routing  metric,  which  aims  for  a  route  with  the 
maximum  mean  detection-probability-to-energy  ratio.  An  algo¬ 
rithm  based  on  combinatorial  optimization  is  also  provided  to 
identify  the  optimal  route.  After  that,  we  consider  the  otergy 
or  detection  i^ormance  constrained  routing  problems,  which 
include  the  minimum  energy  routing  with  guaranteed  detection 
performance  problem  in  Section  IV,  and  the  maximum  detec¬ 
tion-probability  touting  under  constrained  energy  consumption 
problem  in  Section  V.  Linear  integer  programming  formulations 
of  these  problems  are  provided,  together  with  the  methods  and 
algorithms  to  solve  them.  We  provide  a  brief  discussion  in  Sec¬ 
tion  VI,  and  present  the  simulation  results  in  Section  VH.  This 
paper  is  concluded  in  Section  VIII. 

n,  SYSTEM  Description  AND  PROBLEM  Statement 
A,  Sysiem  and  Signal  Model 

Let  us  consider  a  scenario  as  depicted  m  Fig.  1 ,  where  some 
sensors  are  scatteated  over  an  area  to  detect  the  presence  of  an  ob¬ 
ject.  This  scenario  may  correspond  to  monitoring  the  presence 
of  people,  vehicles,  or  military  targets  using  radar-like  sensors* 
that  emanate  specific  eleciromagaetic  signals  into  the  region 
of  interest.  For  this  active  sensing  application,  the  monitored 
space  is  typically  divided  into  many  range  resolution  cells.^ 
l^h  range  cell  could  be  probed  sequentially  in  turn  to  deter¬ 
mine  the  presoic^  of  a  target  by  udng  radar  pulses  that  are  pos¬ 
sibly  launched  by  directional  antamas. 

We  assume  the  position  of  each  sensor  node  is  fixed.^  Bach 
sensor  node  is  assumed  to  be  aware  of  its  own  position  (e.g.,  by 

*Aji  example  of  s«ich  a  rasJa?  sensor  is  ihe  micfc^sowcr  impuiie  radar  [22], 
which  can  be  used  in  wistniciured  sensor  network  environment  wbse  only  lim¬ 
ited  energy  supplies  ane  available.  A  platform  of  tadar-csnabled  sensor  netwodcs 
ia  describe  in  [231 

range  cell  is  the  smalksiii;  range  increment  that  the  radar  k  capable  of 
probing. 

^*T1ie  placement  of  sensors  will  certainly  affect  the  detection  performance  and 
eiieig>*  coiisunipiiorL  Howeva*,  this  is  beyond  scope  of  this  paper,  md  we 
simply  assume  the  sensors  are  already  placed  in  our  ensuing  analysis  and  sim¬ 
ulations.  Interested  readers  can  refer  to  [2^1 
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Fig.  1 .  Routing  in  WSNs  for  signal  detection. 


using  methods  based  on  iriangulation  [25]),  and  has  knowledge 
about  the  location  of  the  sensors  in  its  local  neighborhood.  We 
assume  the  fusion  center  has  information  about  the  geograph¬ 
ical  location  of  each  sensor,  and  will  determine  the  location  to 
be  probed  (for  target  presence)  as  well  as  the  sensors  that  will  be 
involved  in  signal  detection  and  information  transmission.  This 
basically  is  a  destination  initiated  routing,  and  the  routing  algo¬ 
rithm  can  be  implemented  offline,  ahead  of  time,  at  the  fusion 
center.  The  system  will  sequentially  probe  each  location  of  in¬ 
terest,  and  collect  the  information  from  sensor  nodes  lying  on  a 
precomputed  path  in  a  predetermined  order. 

It  is  worth  noting  that  there  exist  several  ways  for  active  sen¬ 
sors  to  illuminate  the  target.  For  example,  in  one  methodology, 
each  sensor  transmits  exactly  the  energy  it  needs  to  obtain  some 
received  signal-to-noise  ratio  (SNR)  if  a  target  is  present  at  a 
certain  location.  This  means  that  the  sensing  unit  should  be  able 
to  dynamically  adjust  its  transmitted  power.  Alternatively,  each 
radar  sensor  could  transmit  the  same  energy  to  the  region  of  in¬ 
terest,  so  that  those  sensors  which  are  closer  to  the  target  will 
have  a  better  observation.  This  might  be  the  situation  when  the 
sensing  unit  is  extremely  simple  and  does  not  allow  power  varia¬ 
tion.  For  brevity,  in  this  paper  we  only  consider  the  second  case, 
but  its  extension  to  the  first  case  is  straightforward.  Therefore, 
let  us  denote  etr  as  the  amount  of  eneigy  with  which  each  sensor 
involved  in  the  signal  detection  uses  to  illuminate  the  target. 
Then  the  energy  received  by  an  antenna  at  a  distance  of  r  me¬ 
ters  away  will  be  ^  •  Ctr  where  ^  is  a  known  constant.  Here  a 
simple  geometric  path  loss  model  [26]  is  assumed,  and  the  path 
loss  is  proportional  to  where  7  is  the  path  loss  exponent, 
an  environment-dependent  constant  typically  between  2  and  4. 

Suppose  the  fusion  center  designates  N  sensors  to  perform 
signal  detection  and  routing.  Thus,  after  matched  filtering  to  the 
known  signal  shape  (which  loses  no  optimality  in  a  Gaussian 


noise  and  interference  environment  [27])  at  sensor  i,  the  obser¬ 
vations  under  the  two  different  hypotheses  are  given  by^ 

y,  =  /\/^  +  ni  (1) 

[ni  Hq 

where  yi  is  the  received  signal,  and  rii  denotes  the  noise  and 
interference.  Moreover,  under  hypothesis  Hit  the  distance  that 
the  active  sensing  signal  traverses,  r,  is  equal  to  where  Ru 
denote  the  distance  between  the  target  (x(f],  2/[<))  and  the  sensor 
(a:(»],y[t]).  i  e., 

Rti  =  v/(xW-x[i])2  +  {yW-y(*])2. 


Since  we  only  focus  on  the  simple  active  sensing  case  in 
this  paper,  for  notational  convenience,  we  define  a  scalar 
0  —  \/^  •  etr/2"^-  Then,  under  //i,  the  received  signal  can  be 
rewritten  as 


Vi  = 


►t/2 


-hni. 


B.  Neyman-Pearson  Delection 


For  the  observation  from  sensor  i,  let  us  define  the  likelihood, 
or  probability  density  function  (pdf)  under  Hi  as  Si{yi)t  and 
the  likelihood  under  Hq  as  /o(s/»)-  Thus,  the  likelihood  ratio  at 
sensor  i  is  simply  given  by 


/■(».) 

/o(w)  fN,{yi) 


where  fsAVi)  denotes  the  pdf  of  the  noise  at  sensor  i.  For  this 
detection  problem,  we  explicitly  make  the  following  assumption 
with  regard  to  the  noise. 

Assumption  I:  The  noise  at  each  sensor  is  assumed  to  be 
statistically  independent,  and  its  probability  distribution  does 
not  depend  on  which  hypothesis  is  true.  □ 

Based  on  this  assumption,  the  joint  probability  of  the  obser¬ 
vations  is  simply  the  product  of  the  individual  probability  den¬ 
sities.  We  define  y  =  [1/1, the  conditional  pdf  of 
the  observation  y  under  Hj  can  be  expressed  as 
N 

i=l 


and  the  log-likelihood  ratio  is  given  by 


Then  under  the  Neyman-Pearson  criterion,  the  optimal  decision 
rule 'd  is  the  solution  to  the  following  constrained  optimization 
problem  [27]: 

max^  PoW 

s.t.  PfW  < 

where  Pni'^)  =  Pr(Choose//i  | //i  true)  and  Pf{^)  = 
Pr{Choose//i  I /fo  true)  are  the  detection  and  false  alarm 

narrowband  signal  model  wiih  random  phase  is  considered  in  [28]. 
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probabilities  when  the  rule  d  is  used,  and  a  is  the  level  of 
the  test.  According  to  the  Neyman-Pearson  lemma  [27]  and 
assuming  no  point  masses  in  the  pdf  of  X(y),  the  optimal 
decision  rule  for  our  hypothesis  testing  problem  can  be  written 
as 

_  /  Decide  iovHu  if  InL(y)  >  v 

\  Decide  (otHq,  if  InL(y)  <  t}  '' 

where  7]  is  uniquely  determined  by  solving  Pf  =  ot. 

It  is  well  known  [27]  that  the  Neyman-Pearson  optimum  test 
(3)  is  a  likelihood  ratio  test.  Since  the  overall  log-likelihood  ratio 
given  in  (2)  is  simply  the  sum  of  the  log-likelihood  ratios  col¬ 
lected  at  each  sensor,  it  follows  that  eaeh  sensor  will  make  a 
partial  contribution  to  the  overall  detection  capability,  which  de¬ 
pends  on  Rii  and  the  noise  properties.  As  the  Gaussian  environ¬ 
ment  is  the  most  important  case  and  is  justified  by  the  central 
limit  theorem,  we  focus  here  on  the  white  Gaussian  noise  case 
by  further  assuming  that  rij  Ar(0,af).  Thus,  the  log-likeli¬ 
hood  ratio  in  (2)  can  be  formulated  into 


lnL(y)  = 


As  a  consequence,  the  likelihood  ratio  test  can  be  written  as 


9  r 


where  we  have  defined  the  test  statistics  g  and  the  new  threshold 
r.  For  a  fixed  set  of  sensors,  the  second  part  of  r  will  be  fixed 
and  known.  The  variable  g  is  actually  a  sufficient  statistic  [20], 
and  when  making  a  decision,  knowing  the  value  of  g  will  be  Just 
as  good  as  knowing  y.  For  notational  convenience,  let  us  define 


N 


i=l 


0^ 


(4) 


where  it  is  noted  that  0^ IcflKl-^  actually  denotes  the  SNR  level 
at  sensor  i.  Then,  invoking  the  model  for  n„  the  hypothesis  pair 
can  be  written  as 


versus 


HQi\g  A/'(0, 0) 


Hi  : 

Thus,  the  false  alarm  probability  is 

Ff  =  Pr(5>r|i/o)  =  l-^>(^) 


(5) 


(6) 


where  is  the  standard  Gaussian  cumulative  distribution 
function,  i.c.. 


Suppose  we  define  the  allowed  level  of  false  alarm  as  Pf  ==  ocy 
then  from  (6),  the  threshold  can  be  computed  by 


r  =  y/ii)  •  \l  —  a). 


Therefore,  the  detection  probability  Po  is  given  by 

FD  =  l-$($-i(l-a)- v^).  (7) 

Since  $(  • )  is  a  monotonically  increasing  function,  we  have  the 
following: 

Proposition  I :  Given  a  specific  value  of  a,  maximizing  Pp 
is  equivalent  to  maximizing  ip,  the  sum  of  SNR  ratios  at  each 
sensor,  or  alternatively:  max  Po  =  max  'ip.  □ 

Hence,  it  is  clear  that  the  term  tj)  essentially  quantifies  the 
detection  performance  gain  earned  across  all  these  N  sen¬ 
sors.  Also,  for  a  fixed  value  of  we  immediately  obtain  the 
following: 

Proposition  2:  Given  a  specific  value  of  \p^  Pu  is  monoton¬ 
ically  increasing  with  respect  to  or.  In  other  words,  in  order  to 
achieve  a  larger  value  of  Pi>,  we  must  also  increase  or.  □ 


C  Energy  Consumption 

In  WSNs,  the  energy  consumption  consists  of  two  parts 
[15]:  one  is  concerned  with  the  communication  between  sensor 
nodes,  and  the  other  with  the  noncommunication  related  energy 
(this  includes  sensing,  signal  processing,  etc.).  We  consider 
these  two  types  of  energy  consumption  as  follows. 

1)  We  assume  that  each  sensor  node  is  able  to  dynamically 

adjust  its  data  transmission  energy  depending  on  the  dis¬ 
tance  over  which  it  transmits  the  signal,  and  that  there  ex¬ 
ists  a  distanee  limit  for  such  information  transmission  (cf. 
Fig.  1 ).  Therefore,  the  transmitting  energy  at  sensor  node  i 
is  where  is  the  distance  between  node  i  and 

node  (i  -I- 1),  which  is  below  the  distance  limit;  ^  is  a  con¬ 
stant;  and  7c  is  the  communication  propagation  coefficient. 
To  maintain  generality,  we  have  assumed  the  propagation 
coefficients  for  data  transmission  and  active  sensing  (e.g., 
7c  and  7)  can  be  different. 

2)  In  Section  II-A,  we  assumed  a  simple  sensing  unit  which 
does  not  allow  power  variation,  and  the  sensing  energy  is 
e^r  However,  it  might  happen  in  large  networks  that  a  node 
participating  in  routing  is  not  able  to  provide  an  essential 
contribution  to  the  overall  detection  performance  when  it 
is  far  away  from  the  target.  This  indicates  that  the  opera¬ 
tions  of  routing  and  sensing  should  be  decoupled  in  prin¬ 
ciple.  Thus,  we  introduce  a  binary  parameter  p,  to  indicate 
whether  node  i  is  involved  in  the  detection  test.  There  are 
several  ways  to  evaluate  this  parameter.  A  simple  method, 
which  we  choose,  is  to  assume  a  sensing  range,  and  let 
only  sensor  nodes  that  are  within  range  of  the  target  lo¬ 
cation  perform  both  sensing  and  routing,  and  the  rest  of 
the  sensors  merely  relay  the  information.  In  this  case,  p,  is 
simply  an  indicator  of  whether  node  i  is  within  the  sensing 
range.5  Therefore,  the  sensing  energy  at  sensor  i  is  pi 

In  addition,  to  avoid  specific  hardware  assumptions,  we 
simply  use  ejp  to  denote  all  other  energy  consumed  at  node 

^It  is  also  possible  to  further  optimize  such  that  it  is  not  tied  to  the  sensing 
range.  This  can  bring  forth  another  optimization  problem  which  however  is  not 
the  focus  of  this  paper. 
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Fig.  2.  Transformed  graph  Q, 


Strongly  connected  network  ^  we  mean  that  for  any  i,j  e 
there  exists  at  least  one  path  which  connects  the  nodes  i  and  j. 
Sinee  targets  are  also  nodes  in  the  graph,  it  is  guaranteed  that 
each  potential  target  location  can  be  reached  by  sensing  signals, 
and  the  information  can  be  transmitted  back  to  the  fusion  center, 
which  is  also  a  node. 

We  denote  the  target  location  (to  be  probed)  as  node  s,  and 
the  fusion  center  as  node  q.  A  path  H  from  s  to  ^  is  a  sequence 
of  N  nonrepeated  nodes  and  connecting  arcs.  Associated  with 
each  arc  (ij)  e  A,  there  are  two  positive  numbers: 

1 )  rpij,  which  specifics  the  detection  performance  gain  earned 
from  the  connection  from  node  i  to  node  j\  or  specifically, 
we  have 


A.. 


t.  This  includes  the  computing  and  signal  processing  en¬ 
ergy  together  with  the  energy  consumed  in  the  transceiver 
circuitry. 

Therefore,  the  overall  energy  consumed  at  the  ith  sensor  can 
be  expressed  as 


2)  Cij  ,  which  denotes  the  energy  consumption  used  by  the 
connection  from  node  i  to  node  j\  or  specifically,  we  have 

^  ^  f  P,  Ctr  +  6Rl°j  +  eip.  i  ^  a 

1  eir  +  e'p  i  =  a. 


fit  —  Pi  Ctr 

We  further  assume  that  there  is  one  sensor  q  placed  at  the  fu¬ 
sion  center,  which  of  course  need  not  perform  data  transmission, 
but  can  operate  at  the  sensing  mode.  Therefore,  the  total  energy 
consumption  for  signal  detection  and  relaying  over  a  route  con¬ 
sisting  of  N  sensors  is 


In  other  words,  the  first  arc  of  each  path  connects  the  target 
to  a  sensor  node,  which  is  a  virtual  path  and  obviously  has  no 
communication-related  energy  consumption. 

Associated  with  any  path  ft  of  the  network  Q  that  connects  s 
to  q,  we  define  two  quantities: 

D{a)^  ^  iPii,  and  EiQ)^  ^ 


N-l 

E  =  (cir  +  c’p)  +  ^  e,. 

1=1 

As  a  special  case,  when  the  probing  location  is  nearby  the  fusion 
center  and  it  suffices  to  detect  the  target  using  only  the  sensor  q 
at  the  fusion  center,  there  will  be  no  actual  routing  involved  and 
the  total  consumed  energy  is  simply  e^r  +  e|p. 

Note  that  what  we  consider  here  is  a  general  model  for  en¬ 
ergy  consumption.  While  it  may  be  desirable  to  consider  a  more 
detailed  energy  consumption  model  in  some  cases,  this  will  in¬ 
evitably  involve  the  consideration  of  hardware  specifications 
that  vary  with  implementation.  To  avoid  such  problems,  we  keep 
our  model  general,  but  we  note  that  our  approach  can  easily  be 
modified  in  an  obvious  way  to  incorporate  alternative  energy 
modeling. 


D.  Graph  Model 

Now  we  introduce  a  graph  model  as  illustrated  in  Fig.  2 
that  will  reflect  the  WSN  topology  (e  g..  Fig.  1).  Moreover, 
this  graph  model  also  incorporates  the  information  on  signal 
detection  performance  and  energy  consumption,  and  will  serve 
as  the  basis  for  our  ensuing  analysis  on  routing  metrics  and 
algorithms. 

Let  Q  =  A)  denote  a  finite  and  strongly  connected  net¬ 
work  defined  by  a  set  M  of  nodes  and  a  set  A  of  directed  arcs. 
The  cardinalities  of  A/'  and  A  are  lAT]  and  |v4|,  respectively.  By 


which  denote  the  detection  gain  and  consumed  energy,  respec¬ 
tively.  For  notational  convenience,  throughout  this  paper,  we  de¬ 
note  5  as  the  (minimum  energy)  shortest  path*^  which  connects 
s  to  9  in  a  graph,  i.e., 

S  =  arg  rnin£(ft). 


in.  Maximum  Path  Detection  Efficiency  Routing 

To  relate  the  detection  performance  to  the  energy  consump¬ 
tion,  we  propose  a  measure  called  the  mean  detection-proba- 
bility-to-energy  ratio,  which  is  defined  as 

A 

E(siy 

According  to  this  measure,  our  objective  is  to  identify  a  path 
which  achieves  the  largest  possible  value  of  A(12),  or  the  max¬ 
imum  path  detection  efficiency.  In  panicular  we  are  only  con¬ 
cerned  with  loopless  paths  or  elementary  paths.  Obviously  this 
metric  aims  at  a  tradeoff  between  the  detection  performance  and 
the  energy  consumption.  To  facilitate  solving  this  problem,  we 

^is  desired  connectivity  can  be  achieved  by  employing  a  topology  control 
technique  in  network  design  (e.g.,  [29]),  which,  however,  is  beyond  the  scope 
of  this  paper. 

^In  the  areas  of  combinatorial  optimization  and  graph  theory,  shortest  path 
usually  refers  to  the  shonest  weighted  path  between  two  nodes  in  an  edge- 
weighted  graph  [30J.  Thus,  in  the  sequel,  we  will  use  minimum-energy  path 
and  shonest  path  interchangeably. 
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choose  to  re-express  our  criterion  as  minimizing 
Thus,  the  problem  can  be  stated  as  follows: 

A\  Given  s,q  £  N,  find  an  optimal  path  which 
satisfies 


...  .  E{Q) 


This  is  a  combinatorial  optimization  problem.  Note  that  when 
.5  =  problem  A  turns  out  to  be  the  minimum  coshto-time 
ratio  cycle  problem,  which  arises  in  an  application  known  as 
the  tramp  steamer  problem  [30,  Ch.  3.  Sec.  13].  This  enables  us 
to  follow  a  similar  parametric  approach  to  that  used  in  [30]  for 
solving  the  minimum  cost-to-time  ratio  cycle  problem.  Such  an 
approach  will  transform  problem  A  into  a  parametric  shortest 
path  problem  [32],  as  described  in  the  following. 

Let  us  denote  p*  as  the  optimal  value  to  problem  A,  i.e., 

♦  A 

p  = 

and  define  p  as  another  nonnegative  scalar.  We  assign  each  arc 
a  modified  length 

hi  =  eij  -  (i^ij 


E{n-) 

£>(«*) 


for  an  arbitrarily  chosen  value  of  p,  and  denote 

9(M,n)=  hi  =  E{i\)  -  nD{iX) 


as  the  length  of  a  specific  path  Also,  let  denote  the  original 
graph  with  new  length  lij  on  each  arc  under  the  parameter  p. 
Adopting  a  similar  geometric  approach  as  in  [31],  we  depict 
q{py  fi)  with  regard  to  p  for  a  given  path  in  Fig,  3,  where 
without  loss  of  generality,  we  assume  5  ^  Since  for  each 
path  ft,  ^{ft)  and  D(ft)  are  strictly  positive,  each  straight  line 
in  Fig,  3  has  a  negative  slope,  and  intersects  the  p  axis  at 


~  D{Q.)  ^ 


Consequently,  the  lower  envelope  of  the  straight  lines,  i.e., 
the  bold  segmented  line  in  Fig.  3  forms  the  shape  of  a  convex 
polygon.  This  bold  boundary  is  the  collection  of  the  shortest 
paths  of  corresponding  to  various  p  values.  Since  the  bold 
boundary  has  a  convex-polygon  shape,  its  intersection  with  the 
p  axis  uniquely  determines  the  optimal  solution  p* . 

Therefore,  algorithms  for  finding  a  solution  to  problem  A  can 
be  designed  to  simply  search  for  the  intersection  of  the  bold 
boundary  and  the  /[x-axis.  Such  searching  procedures  basically 
employ  the  testing  of  successive  trial  values  of  aa  to  find  the 
path  ft^  and  the  optimal  value  p*.  Obviously,  there  are  three 
situations  that  may  exist  with  respect  to  the  path  lengths  q{p^  ft) 
for  the  guessed  value  p  in  the  search  process: 

1)  contains  a  negative  length  path. 

This  case  corresponds  to  q{p.  ft)  <  0,  or  alternatively. 


< 


£;(ft) 

D(ft) 


<  M- 


Thus  /i  is  a  strict  upper  bound  on  p*. 


2)  contains  no  negative  length  path,  but  the  length  of  the 
shortest  path  ft*  is  exactly  zero. 

In  this  case,  the  fact  that  there  is  no  negative  length  path  in 
implies  q{p^  ft)  >  0  or  alternatively, 


M  < 


£(ft) 

Z)(ft) 


for  every  path  ft.  The  faet  that  the  length  of  the  shortest 
path  ft*  is  exactly  zero  implies  that  ^(/z,  ft*)  =  0  or 


p  = 


E{ir) 

D{Q*) 


for  path  ft*.  Thus,  the  guessed  value  p  is  exactly  equal  to 
p* ,  and  ft*  is  the  path  which  gives  the  maximum  detection 
efficiency. 

3)  Qfi  contains  no  negative  length  path,  and  length  of  the 
shortest  path  is  strictly  positive. 

This  case  corresponds  to  ?(//,  ft)  >  0,  or 


p< 


Din) 


< 


for  every  path  ft.  Thus,  the  guessed  value  /x  is  a  strict  lower 
bound  on  p* . 

Based  on  the  preceding  development,  and  adapting  the  bi¬ 
nary  search  procedures  in  [30]  and  [3 1  ],  we  summarize  our  al¬ 
gorithm  for  solving  problem  A  in  Algorithm  1,  where  the  pa¬ 
rameter  e  denotes  the  amount  of  error  that  we  want  to  tolerate 
when  comparing  two  numbers.  Algorithm  1  involves  a  subrou¬ 
tine  for  either  detecting  negative  length  paths,  or  finding  a  non¬ 
negative  length  shortest  path  with  no  repeated  nodes.  An  exact 
implementation  of  such  a  subroutine  has  been  described,  for  ex¬ 
ample  in  [3 1  ],  but  is  very  complicated.  To  lower  the  eomplexity 
of  this  subroutine,  we  propose  an  approximate  approach  for  de¬ 
tecting  negative  length  eycles  instead  of  negative  length  paths, 
whose  effectiveness  is  well  demonstrated  by  our  numerical  re¬ 
sults  in  Section  VII.  To  apply  this  approximate  approach,  we 
can  use  the  Bellman-Ford  algorithm  [33,  Ch.  5]  or  some  other 
improved  negative-cycle  detection  algorithms  (e.g.,  [34]),  for 
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which  the  best  time  complexity  is  0{\A\  •  |A/"|).  Therefore,  the 
time  complexity  of  Algorithm  1  is  approximately  0{\A\  •  jA/'l  • 
log(iiiax(^*j)/c)). 


Algorithm  1:  Maximum  Path  Detection  Efficiency  Routing 

Input:  A  strongly  connected  network  Q  =  {My  A),  and  ejj, 
‘ijiij  for  each  arc  (i,  j)  6  A 

Output:  A  path  il*  of  Q  from  s  loq 
1 :  Use  Dijkstra’s  algorithm  [33,  Ch.  4]  to  obtain  a  shonest 
path  S  from  5  to  ^  of  for  /i  =  0; 

2:  a  ^  E{S)/D{S); 

3:  M  0; 

4:  mark  <—  true; 

5:  While  mark  &Sc  fi  —  fM>  edo 
6:  /i  =  (^  -f-  /2)/2; 

7:  lij  =  Cij  — 

8:  find  a  negative  length  path  O'  (i.e.,  fi')  <  0)  or  a 
nonnegative  length  shonest  path  (i.e.,  O)  >  0). 
from  5  to  ^  of  with  the  current  /x; 

9:  if  any  negative  length  path  ft'  is  found  then 

10:  E{Q')/D{n'); 

1 1 :  else  if  ^(/x,  fi)  =  0  for  the  shortest  path  Cl  then 
12: 

13:  0"  —  Cl; 

14:  mark  false; 

15:  else 

16:  M 

17:  p  ^  m'm{pyE{Cl)/D{Cl)]; 

18:  end  if 

19:  end  while 


IV.  Minimum  Energy  Routing  With  Guaranteed 
Detection  Performance 

The  maximum  path  detection  efficiency  is  meaningful  as  a 
routing  metric,  since  it  reflects  a  certain  relationship  between 
the  detection  performance  and  consumed  energy.  However,  as 
one  drawback,  this  metric  does  not  explicitly  tell  us  what  kind 
of  detection  performance  we  could  achieve  or  how  much  en¬ 
ergy  we  need  to  use.  In  this  section,  instead  of  considering  a 
ratio,  we  are  interested  in  finding  a  route  which  exactly  quanti¬ 
fies  the  detection  performance.  In  other  words,  our  objective  is 
to  identify  a  path  which  minimizes  the  consumed  energy,  while 
maintaining  an  acceptable  level  of  detection  probability. 

To  define  a  desired  level  of  detection  probability,  based  on 
Proposition  /,  it  is  natural  to  employ  a  lower  bound  constraint, 
i.e., 

D(D)  > 

Here  ^  can  be  calculated  from  (7)  for  given  values  of  a  and  Pd 
(e.g.,  a  =  5%,  Pd  =  95%)  as 


±  =  (#-^1  -  a)  -  $-^1  -  Pd)? 
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which  is  clearly  nonnegative.  We  also  impose  an  upper  bound* 
on  given  by  Thus,  this  routing  metric  can  be  described  as 
follows^: 

B:  Given  a^q  €  M,  find  an  optimal  path  il*  which 
satisfies 

=  arg  min  E{il) 

where  11  is  the  set  of  paths  from  s  to  g  for  which  the  detec¬ 
tion  gain  is  bounded  by 

^  <  D{fl)  <  if. 

To  mathematically  formulate  problem  D,  we  adopt  the  idea 
of  network  flow  [33],  and  view  this  problem  as  sending  1  unit 
of  flow  as  cheaply  as  possible  from  node  s  to  ^  subject  to  a  con¬ 
straint  on  the  achieved  performance.  This  interpretation  gives 
rise  to  the  following  integer  programming  formulation,  which 
we  refer  to  as  problem  Bi  : 


Cij  * 

m 

{id)€A 

C  1  i  =  s 

Xji = < 

-1  i  =  9  (9) 

1 0  otherwise 

±< 

(10) 

{id)eA 

Xij  =  {0, 1),  V(i,j)  €  A. 

(11) 

The  first  set  of  constraints  given  in  (9)  are  called  balance 
constraints  which  guarantee  that,  for  all  nodes  on  the  path,  ex¬ 
cept  the  source  and  destination  nodes,  the  inflow  is  equal  to  its 
outflow.  In  the  case  of  the  source  node  {i  =  s),  the  outflow  is 
exactly  one  more  than  its  inflow,  while  in  the  case  of  the  destina¬ 
tion  node  (i  =  q),  its  inflow  is  exactly  one  more  than  its  outflow. 
Equation  (10)  denotes  the  two-sided  weight  constraint,  and  (11) 
is  the  flow  capacity  constraint  of  each  arc.  Each  arc  is  either  in¬ 
cluded  in  the  path  (xij  =  1)  or  not  included  (xij  =  0).  For  con¬ 
venience,  we  denote  /*  as  the  optimal  cost  of  problem  Bi,  and 
define  vectors  x  =  {xij  |  (ij)  €  A),  e  =  {cij  |  (i,  j)  e  A}, 
and  V*  =  {V'ij  I  (»,i)  6  A}. 

Problem  2?i  is  a  constrained  shortest  path  problem,  which  is 
a  special  case  of  a  broader  set  of  optimization  models  known  as 
network  flow  problems  with  side  constraints.  Unfortunately,  this 
problem,  in  the  parlance  of  computational  complexity  theory,  is 
MV-hsrd  [35,  pp.  214],  and  no  polynomial-time  algorithm  is 
available.  In  the  literature,  various  methods  have  been  proposed 
to  solve  this  constrained  shortest  path  problem  (see  a  detailed 
review  in  [37]  or  [38]),  including  Lagrangian  relaxation,  e.g., 
[39]-[43].  This  approach  is  attractive  because  it  can  be  easily 
extended  to  constrained  shortest  path  problems  with  multiple 

^The  upper  bound  ^  appears  redundant  in  this  problem,  but  for  applications 
where  available  resources  are  constrained  (e.g.,  time  or  bandwidth),  this  upper 
bound  would  become  (he  parameter  of  interest.  For  the  purpose  of  completeness, 
we  include  this  upper  bound  here. 

^Ii  is  worth  noting  that  our  formulation  appears  similar  to  constraint-based 
path  selection  in  the  problems  of  qua lity-of- service  (QoS)  routing  uitder  diverse 
network  conditions  or  resource  constraints  (sec  [36]  for  an  overview). 


Authorized  licensed  use  limited  to;  LEHIGH  UNIVERSITY.  Downloaded  on  August  12. 2009  at  19:13  from  IEEE  Xplore.  Restrictions  apply. 


YANG  et  at.:  ENERGY-EFFICIENT  ROUTING  FOR  SIGNAL  DETECTION  IN  WSNs 

side  constrainis,  which  may  arise  when  several  factors  of  in¬ 
terest  are  considered  simultaneously,  e.g.,  delay,  bandwidth  and 
so  on.  Therefore,  in  this  paper,  we  adopt  this  method  to  solve 
the  doubly  constrained  shortest  path  problem  Bi . 

The  Lagrangian  relaxation  approach  basically  relaxes  the 
linear  constraints  by  bringing  them  into  the  objective  function 
with  associated  Lagrange  multipliers  [33).  Let  /ii  >  0  repre¬ 
sent  the  Lagrange  multiplier  for  the  left  inequality  in  (10),  and 
fi2>0be  the  Lagrange  multiplier  for  the  right-hand  inequality 
in  (10).  Then  by  relaxation,  we  obtain  the  resulting  problem 
B2: 


min  ^  [cij  +  ■  Xij  +  -  ti2\p) 


(12) 

(  1  i  =  s 

s.t.  Xij-  Y  ' 

-1  i  =  q 

1 0  otherwise 

(13) 

(14) 

which  is  referred  to  as  the  Lagrangian  subproblem  of  problem 
Bi .  This  program  is  a  minimum  cost  network  flow  problem  of 
sending  1  unit  of  flow  from  sioq  subject  to  the  flow  conserva¬ 
tion  constraints  (13)  with  a  capacity  limit  of  1  on  each  arc  (14). 
This  subproblem  can  be  easily  solved,  although  its  solution  is 
not  necessarily  feasible  for  the  original  problem  Bi.  However,  a 
key  idea  of  Lagrangian  relaxation  is  that  regardless  of  the  choice 
of  pi  and  /i2,  this  subproblem  yields  a  lower  bound  to  the  op¬ 
timal  cost  of  problem  Bi  (i.e.,  /*),  as  shown  by  the  following 
bounding  principle. 

Lemma  1:  Let  X  be  the  set  of  vectors  x  satisfying  constraints 
(13)  and  (14).  For  any  x  6  X,  define  the  Lagrangian  function 


(iJ)eA 

and  let 

-  iTunL(pi,/i2,x). 

(15) 

Then,  it  holds  that 

(16) 

Proof:  See  Appendix  I. 

□ 

It  is  noted  that  the  function  L{pi ,  P2)  defined  in  Lemma  /  is 
concave  and  piecewise  linear,  since  it  is  the  minimum  of  a  finite 
collection  of  linear  functions  of  the  Lagrange  multipliers.  Also 
as  indicated  by  Lemma  7,  to  solve  problem  Bi,  it  is  important 
to  find  the  tightest  possible  lower  bound.  This  involves  finding 
multipliers  and  p2  which  correspond  to 

max  L{hi,(X2) 

S.t.  ^  0* 

This  program  is  referred  to  as  the  Lagrangiofi  dual  problem. 
Solution  to  this  problem  yields  the  sharpest  lower  bound  to  the 
optimal  cost  of  problem  Bi.  Also,  implicated  by  Lemma  /,  we 
immediately  have  the  following  theorem  which  describes  the 
\y^eak  duality. 
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Theorem  1:  The  optimal  solution  of  the  Lagrangian  dual 
problem,  L*  =  max  Z/(/xi,  /i2),  is  always  a  lower  bound  to  the 
optimal  value  of  the  original  problem  i.e.,  L*  <  f*.  □ 

In  other  words,  the  optimal  values  of  pi  and  p2  will  max¬ 
imize  the  Lagrangian  dual  function,  and  give  us  an  optimistic 
estimate  of  the  optimal  value  to  the  original  problem.  Such  an 
estimate  will  be  quite  useful  when  the  main  objective  is  not  to 
find  the  optimum,  but  to  get  a  good  solution.  This  is  another 
attractive  aspect  of  this  Lagrangian  relaxation  approach,  espe¬ 
cially  considering  the  fact  that  the  quality  of  such  an  estimate  is 
guaranteed,  as  shown  by  the  following  theorem: 

Theorem  2:  The  optimal  value  of  the  Lagrangian  dual 
problem,  L*,  is  equal  to  the  optimal  cost  of  the  following  linear 
programming  problem: 

min  e^x 
s.t.  <  ^^x  < 

X  €  conv(X) 

where  conv(X)  denotes  the  convex  hull  of  the  set  X,  i.e.,  the  set 
of  all  convex  combinations  of  the  elements  in  X. 

Proof:  See  Appendix  11.  □ 


Algorithm  2:  Minimum  Energy  Routing  With  Guaranteed 
Detection  Performance 


Input:  A  strongly  connected  network  Q  =  {Jf,  A),  and  Ctj, 
for  each  arc  (i,  j)  6  A 

Output  A  path  Q*  of  Q  from  s  to  q 
1 :  Initialize  the  maximum  lower  bound  found 
Cmax  *  ”001 

2:  Initialize  the  multipliers:  pi  0;  p2  0; 

3:  Initialize  the  iteration  number:  /c  <—  0; 

4:  mark  *—  true; 

5:  while  mark  &  &  k  <  do 
6:  given  the  current  values  of  multipliers  pi  and  use 
shortest  path  algorithms  to  solve  problem  B2,  and  let 
the  solution  be  x; 

7:  compute  the  corresponding  lower  bound:  LB  ♦— 
L(/2i,/22,x); 

8:  \tiUB  50LI?if/c==  0; 

9:  if  LB  >  Craax  then 

10:  Cmax  ^  LB; 

1 1 :  end  if 

12:  if  X  is  a  feasible  solution  to  problem  Bi  then 
13:  UB  ^  e^x; 

14:  end  if 

15:  if  Cmax  =  then 

16:  mark  false; 

17:  else 

18:  calculate  the  subgradients:  ^  ^  h 

+  V^^x; 

19:  compute  the  step  size:  6  =  zo(UB  —  LB)/{g^  -|-  /i^); 
20:  update  the  multipliers  by:  =  max(0,  pi  -|-  tg), 

P2  =  max(0,/i2  +  ^h); 

21:  end  if 

22:  K  K  -f  1; 

23:  end  while 

24:  Resolve  problem  B2  using  p2  that  are  associated 
with  the  resulting  best  lower  bound  Cmaxi 
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To  this  end,  we  summarize  the  algorithm  for  solving  problem 
Bi  based  on  Lagrangian  relaxation  as  shown  in  Algorithm  2, 
which  was  adapted  from  the  subgradieni  method  in  [40],  [41]. 
Here  €  is  a  step  size,  tu  is  a  parameter  satisfying  0  <  <  2, 

and  «niax  is  the  maximum  number  of  iterations  allowed.  This 
algorithm  terminates  when  the  optimal  solution  to  problem  Bi 
is  found,  or  the  iteration  count  k  exceeds  «max*  For  both  cases 
when  the  Bellman-Ford  algorithm  is  used,  the  worst-case  time 
complexity  of  Algorithm  2  will  be  roughly  0(«max  •  Ml '  |*A/^|)- 
However,  in  the  latter  case,  the  best  lower  bound  Cmax  strictly 
less  than  the  optimal  value  of  problem  Bt,  thus  a  duality  gap  is 
said  to  exist,  which  is  also  implied  in  Theorem  I.  Methods  for 
closing  such  a  duality  gap  include  the  path  ranking  approach 
[39],  the  tree  search  procedure  [41],  the  labeling  methods  [42], 
and  so  on.  For  brevity,  we  omit  any  further  discussion  on  this 
issue. 

V.  Maximum  Detection-Probability  Routing  Under 
Constrained  Energy  Consumption 

In  some  applications,  there  might  exist  a  predefined  (or 
estimated)  limited  energy  budget  for  one  single  instance  of  the 
target  detection  and  information  transmission.  Given  such  an 
energy  constraint,  we  may  be  interested  in  finding  a  route  which 
will  achieve  the  maximum  possible  detection  performance,  or 
in  other  words,  make  the  decisions  as  confidently  as  possible. 
Thus,  we  now  consider  the  maximum  detection-probability 
routing  subject  to  a  constraint  on  the  energy  consumption.  This 
is  dual  to  problem  B,  and  can  be  staled  as  follows, 

C:  Given  s^q  €  M,  find  an  optimal  path  ft*  which 
satisfies 


ft*  =  arg  max  B{ft) 

(»j)€n 

where  II  is  the  set  of  paths  from  s  to  ^  for  which  the  energy 
consumption  i.s  bounded  by 

B(ft)  <  £. 


It  follows  that  the  integer  programming  formulation  can  be  ex¬ 
pressed  as  below  (labeled  Ci): 

max  '  ^ij  (12) 


s.t. 


1  i  =  S 

—  I  i  =  q  (18) 

0  otherwise 

(19) 


Xij  =  {0, 1},  y{ij)  e  A.  (20) 


Problem  Ci  turns  out  to  be  a  longest-path  problem  with  con¬ 
straints.  However,  the  longest  path  problem  between  two  nodes 
in  a  weighted  graph  is  already  MV -hardy  since  it  includes  the 
Hamiltonian  path  problem  as  a  special  case,  which  is  one  of  the 
most  well  known  A/^P-hard  problems  [35].  There  exists  some 
work  in  combinatorial  optimization  which  seeks  polynomial¬ 
time  approximation  algorithms  for  longest  paths  (e.g,,  [45]  and 


[46]),  but  these  approaches  are  only  applicable  to  some  specific 
cases.  And,  the  side  constraint  in  (19)  brings  additional  com¬ 
plexity  to  solving  this  longest  path  problem. 

As  it  is  our  inclination  to  formulate  shortest  path  (rather  than 
longest  path)  problems,  for  this  given  problem  Ci  we  propose 
an  algorithm  based  on  a  k-shortest-paih  paradigm  [47],  This 
approach  basically  enumerates,  in  order  of  increasing  length 
(i.e.,  energy  consumption  B(ft)),  all  paths  between  a  and  q  with 
length  smaller  than  the  given  value  E.  Suppose  the  total  number 
of  these  paths  is  M,  We  simply  compare  the  detection  gain  of 
all  these  M  valid  paths,  and  select  the  path  with  the  largest 
value  of  detection  gain.  We  summarize  the  idea  in  Algorithm 
3,  but  it  should  be  noted  that  implementation  of  this  algorithm 
necessitates  a  subroutine  which  can  list  the  shortest  paths  be¬ 
tween  a  and  q  in  order.  An  example  of  such  a  subroutine  is  the 
recursive  enumeration  algorithm  given  in  [48].  If  this  subrou¬ 
tine  is  used,  the  total  time  complexity  of  Algorithm  3  is  about 

0{\A\  +  M■\^/\■log{\A\/m)■ 

This  approach  is  attractive  and  efficient  when  the  network  size 
is  moderate  or  the  difference  between  E  and  E{S)  is  small, 
where  E{S)  is  the  consumed  energy  of  the  shortest  path  S. 
But  in  general,  it  is  hard  to  predict  M,  the  number  of  paths  of 
nondecreasing  length  that  we  need  to  produce  to  satisfy  the  en¬ 
ergy  constraint  in  (19).  Also,  unlike  the  Lagrangian  relaxation 
method,  this  method  does  not  directly  apply  to  the  case  of  mul¬ 
tiple  side  constraints,  and  does  not  provide  a  lower  bound. 


Algorithm  3:  Max-detection-Probability  Routing  With 
Constrained  Energy 


Input  A  strongly  connected  network  Q  =  [M.A)y  and 
il)ij  for  each  arc  {i,j)  E  A 

Output  A  path  ft*  of  G  from  a  io  q 
1:  ({)  *-  0  \ 

2:  i  ^  1; 

3:  mark  *—  true; 

4:  while  mark  do 

5:  find  the  ith  shortest  path  fti  from  a  to  q\nG\ 

6:  if  E(Qi)  <  E  then 
7:  if  B(fti)  >  (f)  then 

8:  ft*  ft*; 

9:  0  ^  Z>(fti); 


10: 

end  if 

11: 

i  <—  z  -f 

12: 

else 

13: 

mark  <— 

14: 

end  if 

15: 

end  while 

VI.  Discussion 

Thus  far,  we  have  focused  on  signal  detection  in  the  Gaussian 
environment.  Unfortunately,  its  extension  to  the  general  non- 
Gaussian  case  is  not  easy,  because  it  requires  the  choice  of  a  spe¬ 
cific  non-Gaussian  pdf,  the  computation  of  the  required  proba¬ 
bilities  (e.g.,  Pf  and  Pd)  is  often  analytically  difficult,  and  it 
seldom  yields  tractable  closed-form  expressions.  As  a  result,  a 
generalization  of  our  problem  formulation  to  the  non-Gaussian 
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environment  is  nontrivial,  although  there  exist  some  very  spe¬ 
cial  cases  (e.g.,  correlation  detection  in  additive  Cauchy  noise 
[27,  p.  84]),  where  it  is  possible  to  obtain  similar  closed-form 
expressions  as  (7)  in  Section  U.  Despite  this,  it  should  be  noted 
that  in  practice,  when  the  pdf  of  the  non-Gaussian  noise  is  un¬ 
known,  the  output  SNR  maximization  seems  natural,  and  can  be 
a  convenient  and  practical  method  for  detection  (e.g.,  [49]).  In 
this  respect,  the  approach  (routing  algorithms)  described  in  this 
paper  could  apply  to  non-Gaussian  environments  as  well,  since 
our  approach  essentially  amounts  to  maximizing  the  sum  of  the 
SNR  ratios  at  each  sensor,  assuming  only  the  second-order  sta¬ 
tistics  of  the  noise  is  known.  Furthermore,  our  approach  is  also 
applicable  to  certain  passive  sensing  scenarios  (cf.  Fig,  1),  for 
example  if  the  system  operates  in  a  similar  way  as  in  the  active 
sensing  case,  by  systematically  detecting  targets  in  geographic 
cells. 

We  can  also  incorporate  node  energy  lifetime,  to  enforce 
route  diversity  and  avoid  unfairly  burdening  some  nodes  (e.g., 
avoiding  hotspot  nodes  [50]).  For  example,  we  can  use  residual 
battery  energy  [51],  as  follows.  Let  bi  denote  the  residual 
battery  capacity  at  node  i,  and  Ci(6i)  be  a  battery  cost  function 
which  denotes  the  reluctance  of  node  i  to  perform  sensing  and 
routing.  Then,  mapping  Ci{bi)  to  the  arc  weight  in  graph  G, 
and  substituting  the  new  weight  for  aj  in  problems  A,  B  and 
C  yields  modified  objective  functions  that  result  in  balancing 
sensor  lifetime  and  the  probability  of  detection.  Consequently, 
all  the  algorithms  provided  in  Sections  UI,  IV,  and  V  would 
still  apply. 


VII.  Simulation  Results 

Now  we  present  simulation  results  for  the  proposed  routing 
metrics  for  detection  in  WSNs.  While,  it  is  highly  desirable  to 
lest  these  algorithms  with  real  data  sets,  due  to  limited  facilities, 
all  of  the  simulation  results  given  in  the  following  are  obtained 
purely  by  specialized  computer  programs  written  in  C-H-  and 
LEDA  [52],  Some  results  based  on  modifying  existing  routing 
protocols  implemented  by  the  network  simulator  ns-2  are  pro¬ 
vided  in  [53]. 

In  our  simulations,  we  consider  WSN  scenarios  where  the 
field  has  a  square  shape  of  size  lO^m  x  lO^m,  50  sensors  are  ran¬ 
domly  distributed  over  this  field,  and  the  fusion  center  is  placed 
at  position  (500,  500),  as  shown  by  a  specific  example  in  Fig.  4. 
We  assume  that  the  fusion  center  wants  to  determine  the  pres¬ 
ence  of  targets  at  different  locations  (e.g.,  stars  in  Fig.  4),  and 
needs  to  find  routes  based  on  our  proposed  metrics  and  algo¬ 
rithms.  We  adopt  7  =  7^  =  2,  ^  ~  1,  and  set  6  =  lO"^, 
(  Ctr  =  0^2'^  =  640  nJ  (nano  joules),  and  Cgp  =  500  nJ  for  all 
i  G  Af.  We  assume  the  data  transmission  range  for  each  sensor 
node  is  250  m,  and  the  sensing  range  is  550  m.  For  simplicity, 
we  assume  (P/a'l  in  (4)  is  the  same  for  every  sensor  node,  and 
further  set  it  to  4  x  10"^,  which  amounts  to  a  6-dB  SNR  ratio  at 
the  distance  of  Rn  =  100  m.  We  use  (equal  to  of  (4)) 
to  quantify  the  detection  performance  over  a  path  Cl.  However, 
for  a  given  false  alarm  probability  Pf  =  a,  the  actual  value  of 
Po  can  be  easily  computed  from  (7).  In  addition,  we  use  our 
propo.sed  approximate  approach  for  step  8  in  Algorithm  1  to 
lower  the  complexity,  and  we  adopt  the  path  ranking  approach 
to  close  the  duality  gap  after  Algorithm  2,  For  the  above  setup, 
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Fig.  4.  Example  of  routing  in  WSNs  for  signal  detection:  (a)  Example  part  A 
and  (b)  Example  pan  B. 


our  computer  programs  can  output  the  results  almost  instantly 
on  an  ordinary  Pentium  computer. 

We  plot  in  Fig.  4  an  example  of  routing  for  detection  in  a  spe¬ 
cific  WSN.  This  example  illustrates  the  different  ways  in  which 
the  routes  from  different  target  locations  to  the  fusion  center  are 
generated  under  our  proposed  routing  metrics  and  some  tradi¬ 
tional  routing  metrics.  In  Fig.  4(a),  we  use  a  to  denote  the  routes 
found  under  the  traditional  minimum-hop  routing  (min-hop);  b 
for  the  routes  under  the  minimum-energy  routing  (min-energy); 
and  c  for  the  paths  under  the  maximum  path  detection  efficiency 
routing  (max-efficiency).  In  Fig.  4(b).  we  denote  d  for  the  routes 
identified  by  the  minimum  energy  rouung  with  guaranteed  de¬ 
tection  performance,  where  we  make  ^  redundant  by  assigning 
it  a  large  value,  and  set  =  2D{S)  (5  is  the  corresponding 
shortest  path  from  a  specific  target  location  to  the  fusion  center). 
Similarly,  we  use  e  in  Fig.  4(b)  for  the  routes  determined  by  the 
maximum  detection-probability  routing  under  constrained  en¬ 
ergy  consumption,  and  set  E  =  2E{S). 

In  Fig,  4(a),  it  is  observed  that  for  some  targets,  the  min-hop 
paths  are  not  the  same  as  the  min-energy  paths,  even  though  both 
of  them  have  the  same  number  of  hops.  This  is  because  min-hop 
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TABLE  1 

Paths  Identified  Under  Different  Routino  Metrics 


;  Type 

j  n  j  £((1)  [  Dfft) 

^0)  Pd 

o 

m»-lwp  ♦JOiiJ  1  1.12 

0J5  '  27.65% 

6 

nw-woiy  ,,  ^3s-4->f  Al9|*.;  |  1.25 

0.30  29  97% 

c 

ttmx-cUKKacy  ->29-«2S~4^v  5.66»iJ  II.JI 

2.00  95.71% 

d 

ikuaksKaasiniiied  -•  13  ->  35  —  4  «  5  4.56 

0.18  1  71.1«% 

€ 

«nefXy<an«tTained  «|->I3-»3U-»3S  — 4~.n— 9  |  1504 

1.66 

paths  may  not  be  unique  and  we  simply  pick  up  the  earliest 
identified  one.  Interestingly  in  this  example,  it  is  seen  that  for 
the  maximum  detection  efficiency  metric,  the  nearest  sensor  to 
the  target  location  is  nearly  always  selected  and  usually  serves 
as  the  actual  source  node  in  the  route.  It  is  intuitively  pleasing 
that  this  route  should  start  with  the  sensor  closest  to  the  target 
location  under  test.  However,  our  joint  optimization  would  not 
necessarily  always  pick  all  the  sensors  that  are  close  to  the  target 
location  [cf.  Fig.  4(a)],  but  instead  will  trade  off  closeness  to 
the  target  location  (i.e.,  those  with  best  SNR)  against  the  energy 
consumption  of  the  route. 

Meanwhile,  to  clearly  illustrate  the  difference  between  those 
routes,  we  consider  the  target  location  (80.  100)  (labeled  as  ^i) 
in  Fig.  4  for  example,  and  list  in  Table  1  all  the  obtained  paths  as 
well  as  the  associated  parameters  of  interest.*®  Also,  by  setting 
a  =  5%  for  instance,  we  compute  the  detection  probability  Pp 
from  (7)  with  respect  to  each  different  path.  It  is  not  surprising 
to  see  that  both  min-hop  and  min-energy  paths  lead  to  inferior 
Po  performance  to  their  counterparts  under  our  proposed  met¬ 
rics,  because  detection  performance  is  not  taken  into  account  by 
either  of  them.  In  addition,  these  results  in  Table  I  also  reveal 
that  our  proposed  three  routing  metrics  are  suitable  to  different 
situations.  To  be  more  specific,  if  one  has  a  stringent  energy 
budget,  he  might  want  to  consider  using  the  maximum  detec¬ 
tion-probability  routing  under  constrained  energy  consumption, 
implemented  by  Algorithm  3.  If  one  cares  more  about  the  de¬ 
tection  performance  (e.g.,  he  wants  to  make  sure  the  decision  is 
at  least  95%  correct),  the  minimum  energy  routing  with  guar¬ 
anteed  detection  performance  implemented  with  Algorithm  2, 
would  be  more  suitable  to  his  need.  These  two  routing  metrics 
are  basically  designed  to  cater  for  different  requirements,  but 
in  essence,  they  are  not  much  different,  and  are  actually  dual 
to  each  other.  For  example,  we  can  set  the  constraints  so  that 
the  two  metrics  provide  the  same  solution.  As  for  the  max -effi¬ 
ciency  routing,  one  might  consider  employing  it  when  there  is 
no  strict  requirement  on  either  enei^y  or  detection  performance, 
but  it  is  desirable  to  squeeze  the  most  performance  return  from 
every  unit  of  energy  put  in  both  sensing  and  routing. 

Funher,  to  illustrate  the  affect  of  the  lower  bound  rp  over  the 
route  selection  when  solving  problem  J3i,  we  condu^ a  simple 
experiment  for  the  target  location  (920,  920)  (referred  to  as  S2) 
at  the  upper  right  comer  of  Fig.  4.  In  this  instance,  we  assign 
^  different  values,  i.e.,  jp  =  D{S),  1.5D(5),  2D{S),  2,5D{S), 
and  9. 5/7(5).  For  each  value  of  rp,  we  resolve  problem  Bi ,  and 
record  the  selected  path.  The  results  are  shown  in  Table  II,  in 
which  we  also  provide  the  values  of  Pq  for  a  =  5%.  Intuitively, 
different  lower  bounds  tp  may  lead  to  different  paths.  For  ex¬ 
ample,  when  ^  =  Z)(5),  the  constrained  shoncst  path  is  simply 
the  min-eneigy  path;  when  ^  =  2.5/7(S).  the  resulting  path, 
for  this  example,  happens  to  be  the  max -efficiency  path.  How¬ 
ever,  it  should  be  noted  that  due  to  the  discrete  characteristics  of 

*®Thc  unit  for  ihc  detecuon  efficiency  A(D)  is  l/fiJ. 


TABLE  11 

Paths  Found  by  Varying  ^  in  Problem  Dx 
FOR  Target  Location  (92^r,  920);  a  =  5% 


n 

m) 

D(n) 

Pd 

l.O 

*  $2  7  -*2  —*  q 

4.39/iJ 

1.28 

30.38% 

1.5 

JS2  — *  43  -♦  2  — ►  g 

4.44  fiJ 

2.27 

44.50% 

2.0 

92  — 9  —  38  —  27  — 

5.60 /ij 

2.80 

51.14% 

2.5 

•a  —  41  —  43  -  2  — 

5.67  fiJ 

11.62 

96.12% 

9.5 

52  —  41  — 9  —  28  —  2— 

6.87 

12.28 

96.85% 

•  be 
Knrijo'  OiiiMun|iUMi  A’llIJ  (iu  /il) 


.Mi'iiii  /*/,  (ill  *X)  .Mphji  IVirriifia  Effirifiiry  A(l| 

Fig.  5.  Average  parameters  of  interest. 


problem  Z?i ,  a  change  to  ^  does  not  necessarily  lead  to  a  change 
in  the  resulting  path.  For  instance,  when  we  set  ^  =  9/7(5),  the 
output  of  the  algorithm  remains  to  be  the  max-cfficiency  path  for 
the  case  in  Table  11.  A  similar  phenomenon  also  occurs  when  we 
solve  Cl  under  different  energy  constraints  (i.e.,  E).  We  omit 
further  results  and  discussion  for  brevity. 

We  also  investigate  the  statistical  performance  of  the  three 
routing  metrics:  min-hop,  min-energy,  and  max-efficiency, 
which  are  basically  parameter  free  and  do  not  depend  on 
any  side  constraints.  In  our  experiments,  we  generate  25  090 
random  sensor  placements,  and  one  random  target  location  for 
each  sensor  placement  is  tested.  We  identify  the  corresponding 
routes  and  calculate  the  performance  metrics  of  interest.  In 
particular,  we  compute  Pq  given  a  =  5%,  and  use  it  instead 
of  D{fl)  to  represent  the  detection  performance.  We  average 
these  results,  and  plot  them  in  Fig.  5.  For  a  fair  comparison,  we 
also  plot  the  empirical  cumulative  distribution  function  (CDF) 
of  the  energy  consumption  detection  probability  Pp, 

and  detection  efficiency  X{il)  in  Fig.  6.  From  both  figures, 
it  is  shown  that  our  max-efficiency  routing  uses  about  1  ^ij 
more  energy  than  the  other  two.  This  is  not  surprising  because 
min-energy  routing  consumes  the  least  eneigy  and  min-hop 
routing  closely  resembles  it  in  our  experiment  set-up.  However, 
the  superiority  of  the  proposed  max -efficiency  metric  can  be 
easily  verified  from  Figs.  5  and  6  by  examining  the  detection 
probability  and  efficiency.  For  example,  as  shown  in  Fig.  6(b), 
almost  50%  of  the  resulting  paths  under  our  metric  give  a 
detection  probability  Pp  which  is  larger  or  equal  to  90%, 
while  only  around  3%  of  the  paths  under  either  min-hop  or 
min-energy  metric  have  such  performance.  Therefore,  using 
a  path  which  consumes  a  bit  more  energy  than  a  path  under 
the  traditional  methods,  we  actually  have  utilized  every  unit 
of  consumed  energy  much  more  efficiently.  This  demonstrates 
achieving  the  desired  tradeoff  between  energy  consumption 
and  detection  performance. 
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(a) 


(c) 

Fig.  6.  Empirical  CDF  for  paramcicrs  of  imeresl.  (a)  Empirical  CDF  of  energy 
consumption  £((7)  (in/j  J),  (b)  Empirical  CDF  of  Pp  in  %,  (c)  Empirical  CDF 
of  detection  efficiency  A(n). 


Vlll.  Conclusion  and  Future  Work 

For  a  target  location  problem  of  great  interest,  we  formulated 
the  problem  of  energy-efficient  routing  for  detection  in  a  wire¬ 
less  sensor  network  under  the  Neyman-Pearson  detection  crite¬ 
rion,  which  relates  both  the  energy  consumption  and  detection 
performance  (in  terms  of  the  probabilities  of  detection  and  false 
alarm)  in  routing.  We  also  proposed  three  routing  metrics  which 
strive  for  a  balance  between  the  consumed  energy  and  detec¬ 
tion  performance.  With  respect  to  each  different  routing  metric, 
we  formulated  the  corresponding  routing  problem  into  a  com¬ 
binatorial  optimization  program,  and  further  provided  an  effi¬ 
cient  algorithm  for  finding  the  optimum  solution  based  on  tech¬ 
niques  from  operations  research.  A  specific  simulation  example 
was  presented,  which  illustrates  the  manner  in  which  the  paths 
are  selected  under  these  proposed  different  routing  metrics.  We 
also  provided  an  example  which  indicates  the  impact  of  the 
side  constraint  on  route  selection,  and  further  demonstrates  the 


flexibility  of  our  algorithms  in  balancing  various  factors.  More¬ 
over,  Monte  Carlo  simulation  results  based  on  generating  a  large 
number  of  random  topologies  were  also  given,  demonstrating 
how  the  maximum  path  detection  efficiency  metric  achieves  su¬ 
perior  detection  performance  and  energy  consumption  tradeoffs 
over  the  traditional  min-hop  and  min-energy  routing  metrics. 

Our  work,  however,  can  be  further  improved  and  extended. 
For  example,  one  could  easily  imagine  using  several  parallel 
routes  for  robustness.  Further,  the  detection/routing  idea  could 
be  generalized  to  include  target  tracking.  It  might  also  be  in¬ 
teresting  to  consider  an  extensive  study  of  the  network  lifetime 
jointly  with  the  detection  performance  in  the  routing  problem. 
This  in  turn  will  require  more  statistical  data  with  regard  to  the 
flow  status  of  each  node,  which  we  did  not  assume  in  this  paper. 
Further,  it  would  be  desirable  to  develop  distributed  routing  al¬ 
gorithms  for  multitarget  detection  and  tracking. 


Appendix  1 
Proof  of  Lemma  1 

Proof:  Using  vectors  x,  e,  and  the  Lagrangian  function 
can  be  reformulated  into 

=  e^x  +  (/i2  -  +  (/*1^  - 

and  the  two-side  weight  constraint  in  (10)  can  be  written  as  ^  < 
Thus,  we  have  the  following: 

^(M1,M2) 

=  min{e^x  +  {/i2  -  +  (/iiV'  -  M2^)} 

x€X  — 

<  min  {e^x 

+  +  M2(V’^x  -  i))} 

<  min  e^x 

x6X,^<^^x<iii« 

where  the  equality  follows  from  the  definition  (15);  the  first  in¬ 
equality  follows  because  the  minimum  of  the  right-hand  side  of 
this  inequality  is  taken  over  a  subset  of  X  (x  has  to  satisfy  the 
doubly  constraint  as  well);  and  the  second  inequality  follows 
because  of  the  nonnegativity  of  /ii  and  fi2.  Since  the  right-hand 
side  of  the  second  inequality  denotes  the  optimal  value  of  the 
problem  Di,  i.e., 

/*  =  min^  e^x, 

x  6 ,  X ,  ^  X  <-0 

we  immediately  obtain  (16).  □ 


Appendix  11 
Proof  of  Theorem  2 

Proof:  By  the  definition  in  (15),  £(/2i,/i2)  is  linear  in  x, 
thus  the  optimal  cost  remains  the  same  if  we  allow  convex  com¬ 
binations  of  the  elements  of  X.  As  a  result,  we  have 

i(Mi,M2)=  min  L(Mi,/i2.x) 
x€conv(X) 
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and  further 


L*  =  max  min 

x€coav(X) 

Since  the  set  X  has  finite  elements,  conv(X)  is  a  bounded  poly¬ 
hedron.  Let  X*,  fe  6  /C,  be  the  extreme  points  of  conv(X),  where 
/C  =  {1, . . . ,  AT}  is  a  finite  set  of  integers.  Then,  for  any  fixed 
values  of  iii  and  ^2*  we  have 

L{ni,fi2)  =  nuuL(/xi,M2,x*) 

which  follows  from  the  faei  that  for  a  linear  programming 
problem  of  minimization,  there  exists  an  extreme  point  which 
is  optimal  if  the  optimal  cost  is  not  unbounded  [44,  Theorem 
2.8]. 

Therefore,  the  Lagrangian  dual  problem  is  equivalent  to  and 
has  the  same  optimal  value  L*  as  the  problem 

max  minL(/ii,^2jX*^) 

S.t. 

By  introducing  an  auxiliary  variable  z,  the  above  problem  can 
be  equivalently  reformulated  into  the  following: 

max  z 

s.t.  z  -h  +  /^2(^  - 

<  €  AC 

>  0. 


Taking  the  linear  programming  dual  of  the  above  problem,  and 
applying  strong  duality  [44,  Theorem  4.4],  we  can  see  that  L* 
is  equal  to  the  optimal  cost  of  the  following  problem: 


Ofjt  >  0,  A:  G  AC. 


Since  conv(X)  denotes  the  set  of  all  convex  combinations  of  the 
elements  in  X,  we  have 


coav(X) = 


a/t  =  1,  ajt  >  0,  A:  G  AC 

kefc 


Therefore,  the  foregoing  linear  programming  dual  can  be  refor¬ 
mulated  into 

min  e^x 
s.t.  jp  <  %lP‘x  <  \p 
X  G  conv(X) 


whose  optimal  cost  is  equal  to  L* . 


□ 
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